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Abstract 


To identify host factors relevant for SARS-coronavirus (SARS-CoV) replication, we performed an siRNA 
library screen targeting the human kinome. Protein kinases are key regulators of many cellular functions 
and the systematic knockdown of their expression should provide a broad perspective on factors and 
pathways promoting or antagonizing coronavirus replication. In addition to 40 proteins that promote 
SARS-CoV replication, our study identified 90 factors exhibiting an antiviral effect. Pathway analysis 
grouped subsets of these factors in specific cellular processes, including the innate immune response and 
the metabolism of complex lipids, which thus appear to play a role in SARS-CoV infection. Several 
factors were selected for in-depth validation in follow-up experiments. In cells depleted for the B2 subunit 
of the coatomer protein complex (COPB2), the strongest proviral hit, we observed reduced SARS-CoV 
protein expression and a >2-log reduction in virus yield. Knockdown of the COPB2-related proteins 
COPB1 and Golgi-specific brefeldin A-resistance guanine nucleotide exchange factor | (GBF1) also 
suggested that COPI-coated vesicles and/or the early secretory pathway are important for SARS-CoV 
replication. Depletion of the antiviral double-stranded RNA-activated protein kinase (PKR) enhanced 
virus replication in the primary screen, and validation experiments confirmed increased SARS-CoV 
protein expression and virus production upon PKR depletion. In addition, cyclin dependent kinase 6 
(CDK6) was identified as a novel antiviral host factor in SARS-CoV replication. The inventory of pro- 
and antiviral host factors and pathways described here substantiates and expands our understanding of 


SARS-CoV replication and may contribute to the identification of novel targets for antiviral therapy. 
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Importance 


Replication of all viruses including SARS-coronavirus (SARS-CoV) depends on and is influenced by 
cellular pathways. Although substantial progress has been made in dissecting the coronavirus replicative 
cycle, our understanding of the host factors that stimulate (proviral factors) or restrict (antiviral factors) 
infection remains far from complete. To study the role of host proteins in SARS-CoV infection, we set out 
to systematically identify kinase-regulated processes that influence virus replication. Protein kinases are 
key regulators in signal transduction, control a wide variety of cellular processes, and many of them are 
targets of approved drugs and other compounds. Our screen identified a variety of hits and will form the 
basis for more detailed follow-up studies that should contribute to a better understanding of SARS-CoV 
replication and coronavirus-host interactions in general. The identified factors could be interesting targets 
for the development of host-directed antiviral therapy to treat infections with SARS-CoV or other 


pathogenic coronaviruses. 
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Introduction 


Positive-stranded RNA (+RNA) viruses interact with the infected host cell at many levels during their 
replicative cycle, and thus far numerous host cell proteins that influence virus infection have been 
identified (1, 2). These include, for example, host factors recruited by the virus during the various stages 
of its replicative cycle and those involved in the host’s defense against virus infection. Such proteins may 
constitute interesting targets for the development of novel antiviral strategies, as drug resistance is less 
likely to develop when cellular rather than viral functions are targeted. Antiviral drug resistance is a 
serious problem, in particular when combating RNA viruses, due to their high mutation rate and potential 
for rapid adaptation. 

Systems biology approaches have been instrumental in advancing our knowledge of the proteins 
and cellular pathways that influence +RNA virus infection. For example, systematic functional genomics 
screens using small interfering RNA (siRNA) libraries have identified numerous host proteins with a role 
in the replication of important human pathogens like West Nile virus (3), Dengue virus (4, 5), human 
immunodeficiency virus | (6), hepatitis C virus (7-12), and influenza virus (8, 13, 14). For coronaviruses a 
number of relevant host proteins have previously been described ((15-17), and reviewed in (2, 18)), but 
the use of siRNA screens to systematically identify such factors has not been reported thus far. 

Coronaviruses, and some other members of the order Nidovirales (19), have the largest RNA 
genomes known to date (25-34 kb (20)) and the complexity of their molecular biology clearly 
distinguishes them from other +RNA virus groups. Although infection with most established human 
coronaviruses is associated with relatively mild respiratory symptoms (21, 22), the 2003 outbreak of 
severe acute respiratory syndrome (SARS) highlighted the potential of zoonotic coronaviruses to cause 
lethal disease in humans. The emergence of SARS-coronavirus (SARS-CoV), which likely originated 
from bats, initiated an outbreak that affected about 8,000 humans, with a mortality rate of approximately 


10% (23). Strikingly, a similar outbreak of coronavirus-induced severe respiratory disease has been 
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developing in a number of Arab countries since April 2012, with ~420 of the >1100 confirmed cases 
having a fatal outcome by April 2015 (http://www.who.int/). The causative agent, Middle East Respiratory 
Syndrome-coronavirus (MERS-CoV), was identified as a previously unknown member of the 
betacoronavirus subgroup 2c (24, 25). These recent developments stress the importance of developing 
antiviral approaches to combat coronavirus infections and highlight the relevance of the systematic 
dissection of coronavirus-host interactions. 

SARS-CoV RNA synthesis, like that of many +RNA viruses (26), takes place at virus-induced 
membrane structures (27, 28), which in this case comprise a reticulovesicular network (RVN) of modified 
endoplasmic reticulum ((28) and reviewed in (29)). The viral replication and transcription complexes 
(RTCs) are associated with this RVN, which is thought to create a suitable microenvironment for RNA 
synthesis and possibly also provides protection against cellular antiviral activities. The biogenesis of the 
RVN, and the functional details of the RTC, in particular the role of cellular factors and pathways, are far 
from understood. 

Previous studies addressed coronavirus-induced immune responses, as well as a number of 
specific interactions between coronaviruses and the antiviral immune response (reviewed in (2)). Several 
immune evasion mechanisms were attributed to protein functions that are either conserved across CoVs or 
specific for certain CoV lineages. Proteins such as non-structural protein 1 (nsp1; (30), the nsp3 papain- 
like proteinase (31), the nsp16 2’-O-methyltransferase (32), the nucleocapsid (N) protein (33), and the 
products of SARS-CoV ORFs 3b and 6 (34-37) have been reported to interfere with interferon (IFN) 
induction and/or signalling. In addition, the SARS-CoV E protein has been shown to manipulate the 
cellular stress response in cell culture, including the unfolded protein response and apoptosis (38). 

To gain more insight into the role of host factors in the SARS-CoV replicative cycle, we set out to 
systematically identify kinase-regulated cellular processes that influence virus replication. Protein kinases 
are key regulators in signal transduction and control a wide variety of cellular processes. Thus, assessing 


their relevance for virus replication can provide a broad perspective on factors and pathways relevant for 
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SARS-CoV replication, as illustrated by previous studies identifying cellular kinases as host factors 
influencing various stages of the replicative cycle of other +RNA viruses (5, 10, 11, 39, 40). 

In this study, we have screened an siRNA library that targets the cellular kinome (779 genes) and 
identified 40 proviral and 90 antiviral factors whose depletion significantly reduced or enhanced SARS- 
CoV replication, respectively. Pathway analysis grouped several subsets of hits in specific cellular 
pathways, suggesting that these play an important role in the SARS-CoV-infected cell. Two strong hits 
from the siRNA screen, the proviral B2 subunit of the coatomer complex (COPB2) and the antiviral 
double-stranded RNA-activated protein kinase (PKR), were selected for independent validation and 
follow-up analysis, which confirmed their importance for SARS-CoV replication. In addition, several 
other hits from the primary screen were evaluated, and the relevance of the antiviral factor CDK6 and the 
proviral factor PRKC1 could be confirmed. Our data offer a glimpse into the complex interplay between 
SARS-CoV and its host cell, and provide a basis for in-depth studies that will enhance our understanding 


of coronavirus replication and coronavirus-host interactions. 


Materials and methods 


Cell culture, compound, viruses, and virus titration — 293/ACE2 (41) and Vero E6 cells were cultured 
as described previously (42). Although 293/ACE2 cells have been described as a human 293 cell-derived 
cell line (41), our recent work established that these cells actually must have originated from a non-human 
primate species that is closely related to the rhesus monkeys Macaca mulatta and Papio Anubis (43). Cells 
were infected with SARS-CoV strain Frankfurt-1 (44) or GFP-expressing recombinant SARS-CoV 
(Urbani strain) (45) as described previously (42). Sodium aurothiomalate (ATM; Sigma cat. nr. 157201) 
was dissolved in PBS and stored as 100 mM stock at -20°C. Virus titrations were performed essentially as 
described before (46). All work with infectious wild-type (wt) SARS-CoV and SARS-CoV-GFP was 
performed inside biosafety cabinets in a biosafety level 3 facility at Leiden University Medical Center. 
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siRNA library and transfection reagents — The ON-TARGETplus SMARTpool Protein Kinases siRNA 
Library that targets the mRNAs of 779 genes, comprising the complete human kinome and some 
additional targets, was obtained from Dharmacon. Each individual siRNA SMARTpool consisted of four 
siRNAs targeting the same gene. A non-targeting (scrambled) siRNA (cat. nr. D-001810-10; Dharmacon) 
served as a negative control and a GAPDH-targeting siRNA (cat. nr. D-001830-10; Dharmacon) was used 
to routinely monitor transfection and knockdown efficiency. Stock solutions (2 M) of siRNA 
SMARTpools were prepared by dissolving 0.5 nmol of an siRNA SMARTpool in 250 pl of Ix siRNA 
buffer (Dharmacon), according to the manufacturer’s instructions. Using a 96-well pipettor (Rainin 
Liquidator 96), the contents of the siRNA library master plates was aliquoted into volumes appropriate for 
individual screening experiments. The resulting sets of ten deep-well 96-well library plates (Greiner Bio- 


One) were stored at -80°C until further use. 


siRNA library screening and validation - In each siRNA screen, 293/ACE2 cells in 96-well plates 
containing ~10* cells per well were transfected with a 100-11 mixture containing 100 nM siRNA, 0.2 pg 
DharmaFECT1 (Dharmacon), OptiMEM (Invitrogen), and antibiotic-free cell culture medium, 
supplemented with 8% fetal calf serum (FCS) and 2.5 mM L-Glutamine, according to Dharmacon's 
instructions. Transfection mixes were prepared in the ten deep-well 96-well plates that together contained 
the complete library of 779 siRNA SMARTpools (see above). Using the contents of these library plates, 
we transfected 293/ACE2 cells in black (3 wells per target) and transparent 96-well plates (3 wells per 
target). For a schematic representation of the experimental set-up, see Fig. 2. Transfection of individual 
siRNAs (ON-TARGETplus siRNAs; Dharmacon) targeting CDK6 (cat. nr. LU-003240-00), MAP2K1 
(cat. nr. LU-003571-00), MAP2K3, (cat. nr. LU-003509-00), PKR (cat. nr. LU-003527-00), or siRNA 
SMARTpools targeting COPB1 (cat. nr. L-017940-01) and GBF1 (cat. nr. L-019783-00) was performed 


as described previously (42). Twenty-four hours post transfection (p.t.), the medium was replaced, and 
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cells were incubated for another 24 h at 37°C. At 48 h p.t., cells were infected with SARS-CoV-GFP at an 
MOI of 10, and 24 h later they were fixed with 3% paraformaldehyde (PFA) in PBS. GFP expression was 
quantified by measuring fluorescence in a 96-well plate reader (Berthold Mithras LB 940), using 
excitation and emission wavelengths of 485 and 535 nm, respectively. The fluorescence in wells 


containing mock-infected cells was used to correct for background signal. 


GAPDH and cell viability assays - At 48 h p.t., GAPDH enzyme activity in lysates of siRNA-transfected 
cells was measured using the KDalert™ GAPDH Assay Kit (Ambion) according to the manufacturer's 
instructions. Possible cytotoxic effects of siRNA transfection were analyzed (in triplicate) at 48 h p.t., 
using the CellTiter 96® AQueous Non-Radioactive Cell Proliferation Assay (Promega). After 90 min, the 
reaction was terminated by the addition of 25 ul of 10% SDS and absorbance at 490 nm (A4oo) was 


measured using a 96-well plate reader (Berthold). 


Data analysis — Raw data from GFP fluorescence and cell viability measurements were analyzed per 
individual screen with the Bioconductor/R package CellHTS2 (47) with minor modifications (see Results 
section for details). Average GFP expression (n=3) and cell viability were calculated and normalized to 
the signals of scrambled siRNA-transfected (control) cells. A two-sided one-sample Student’s ¢ test was 
used on the log»-transformed normalized values to determine the significance (p < 0.05) of the changes in 
GFP expression caused by siRNA transfection. The siRNA transfection was considered non-cytotoxic 
when the normalized cell viability assay readings (Ao) were above 0.85 (p < 0.05). Significance was 
determined using a one-sided one-sample Student’s ¢ test on the log»-transformed normalized values using 


p< 0.85 as the null hypothesis. 


Gene silencing using lentivirus-expressed shRNAs — Vectors for expression of short hairpin RNAs 


(shRNAs) targeting human COPB2 (cat. nr. TRCN-065114; accession nr. NM_004766) or expression of a 
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non-targeting (scrambled) control shRNA (cat. nr. SHC-002) were picked from the MISSION TRC-1 
library of shRNA-expressing lentiviruses (Sigma) and lentivirus stocks were prepared according to the 
manufacturer’s instructions. Lentivirus particle titers were determined using a p24 ELISA (Zeptometrix) 
according to the manufacturer's instructions. Wells (4 cm’) containing 8 x 10* 293/ACE2 cells were 
transduced with shRNA-expressing lentiviruses at an MOI of 3 in culture medium containing 8 pg/ml 
polybrene, and after 24 h fresh medium was given. At 72 h p.t., cells were infected with wt SARS-CoV or 


SARS-CoV-GFP (MOI 0.01), and depletion of COPB2 was validated by Western blotting. 


Protein analysis and antibodies — Total cell lysates were prepared in 4x Laemmli sample buffer (100 
mM Tris-HCl, pH 6.8, 40% glycerol, 8% sodium dodecyl sulfate (SDS), 40 mM DTT, 0.04 mg/ml 
bromophenol blue), after which samples were heated at 95°C for 15 min. Following SDS-PAGE, proteins 
were transferred to Hybond-LFP membranes (GE Healthcare) by semi-dry blotting, and membranes were 
blocked with 1% casein in PBS containing 0.1% Tween-20 (PBST). The following antisera against 
cellular proteins were used: rabbit anti-PKR (cat. nr. 610764; BD Biosciences), goat anti-COPB2 (sc- 
13332; Santa-Cruz), rabbit anti-CDK6 (sc-177; Santa Cruz), rabbit anti-MAP2K1 (710446; Life 
Technlogies), rabbit anti-MAP2K3 (sc-961; Santa Cruz), and mouse monoclonal antibodies against B- 
actin (A5316; Sigma) and the transferrin receptor (TfR; cat. nr. 13-6890; Invitrogen). Rabbit antisera 
against SARS-CoV nsp8 and N protein (28, 48) were used to analyze viral protein expression. After 
overnight incubation with the primary antibody, membranes were probed with biotinylated secondary 
antibodies (rabbit anti-goat, swine anti-rabbit, or goat anti-mouse) for 1 h at RT, after which a tertiary 
mouse anti-biotin-Cy3 antibody was used to visualize protein bands using a Typhoon 9410 scanner (GE 


Healthcare). 


Canonical pathway analysis - The Ingenuity Pathway Analysis IPA™) package was used to assign hits 


to canonical cellular pathways. The significance of the association between the dataset and the respective 
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pathways was determined in two ways: (i) the number of molecules from the dataset that mapped to a 
specific pathway divided by the total number of molecules in that canonical pathway (the higher the 
percentage of hits identified in a specific pathway, the higher the likelihood it plays a role in the viral 
replicative cycle); (ii) Fisher’s exact test was used to determine the probability that the association 


between the genes in the dataset and the canonical pathway is explained by chance alone. 
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Results 


Development of an siRNA screening protocol for host factors involved in nidovirus replication. A 
robust protocol was developed to assess the effect of systematic knockdown of individual host kinases on 
the replicative cycle of SARS-CoV (this study) and the distantly related arterivirus equine arteritis virus 
(EAV; K. F. Wannee, A. H. de Wilde et a/., manuscript in preparation), which was applied to the 
screening of a commercial human kinome-directed siRNA library (779 targets). We performed our siRNA 
screens in 293/ACE2 cells (41), which express the SARS-CoV receptor angiotensin-converting enzyme 2 
and, in contrast to other cell lines tested, were found to be permissive to a combination of siRNA 
transfection and infection with either SARS-CoV or EAV. This property facilitated direct comparative 
studies between these two distantly related nidoviruses. Unfortunately, after completion of the siRNA 
screens, it was discovered that these cells are not of human origin, but have most likely originated from an 
Old World monkey closely related to Papio Anubis and Macaca mulatta (43). Nevertheless, because the 
sequence identity between the human genome and that of several Old World monkeys is 94% (49) and 
because pools of four siRNAs were used for each target, we believe that the consequences of the 
misidentification of this cell line are limited, although the chance of false-negative hits may have been 
somewhat increased (49). Infection of 293/ACE2 cells with SARS-CoV-GFP at an MOI of 10 yielded a 
robust and readily detectable GFP signal at 24 h p.i. (Fig. 1A). The GFP signal was stronger at 28 and 30 h 
p.i., indicating that it had not yet reached a plateau at 24 h p.i. (Fig. 1A). We therefore chose the latter 
time point to fix cells and measure GFP fluorescence, as it should also allow the identification of antiviral 
factors whose knockdown would increase reporter gene expression. The 293/ACE2 cells could be 
efficiently transfected with siRNAs, as illustrated by a consistent ~75% reduction of GAPDH activity at 
48 h p.t. using an siRNA SMARTpool targeting the GAPDH mRNA (Fig. 1B; white bars). No change in 
cell viability was detected by 48 h p.t. following transfection with either a scrambled siRNA or the 


GAPDH-specific siRNA (Fig. 1B; grey bars). When these cells were subsequently infected with SARS- 
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CoV-GFP (MOI 10), no significant differences in GFP expression were observed at 24 h p.i. compared to 
control cells that had not been transfected with siRNAs. This demonstrated that the siRNA transfection 


procedure per se did not adversely affect SARS-CoV-GFP replication (Fig. 1B; black bars). 


Kinome-wide siRNA screens for host factors involved in SARS-CoV replication. A human kinome- 
directed siRNA screen was performed to identify host cell kinases that affect SARS-CoV-GFP replication, 
according to the experimental set-up outlined in Fig. 2. For each independent siRNA screening 
experiment, we used a set of ten 96-well library plates, each containing approximately 80 specific siRNA 
SMARTpools and several controls. Transfection mixes (final concentration of 100 nM siRNA) were 
prepared in these library plates and their contents was used to transfect - per library plate - 293/ACE2 cells 
in three black and three transparent 96-well plates. Forty-eight hours after siRNA transfection, the black 
plates were infected with the SARS-CoV-GFP reporter virus (MOI 10), and at 24 h p.i. GFP expression 
was measured by fluorometry. At the moment of infection, the transparent plates were used to monitor 
(potential) cytotoxic effects of siRNA transfection using a colorimetric cell viability assay. The complete 
siRNA screen, i.e. the viability controls (in triplicate for each siRNA SMARTpool) and the quantitation of 
SARS-CoV-driven GFP expression (in triplicate), was repeated in three independent experiments. The 
data, obtained from a 96-well plate reader, was processed with the Bioconductor/R package CellHTS2 as 
described (47). Experimental controls were assigned, and the NPI method (normalized percent of 
inhibition) was used to normalize GFP fluorescence values to those of scrambled siRNA-transfected cells, 
and to correct for plate-to-plate variation. Subsequently, the GFP data were transformed to a multiplicative 
scale (the value obtained using scrambled siRNA-transfected cells was set to 1). Next, the results for each 
replicate library screen were summarized and used for further data analysis, including the assignment of 
GenelIDs to each well. Finally, the data of the three independent library screens were combined and 


summarized. 
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Host cell kinases were considered to have a proviral effect when their siRNA-mediated 
knockdown reduced the GFP signal (negative score values) and kinases were considered antiviral when 
the GFP signal increased upon their knockdown (positive score values). Graphical representations of the 
hit distribution per plate were visually inspected in order to minimize the chance of false positive or false 
negative hits due to major (technical) artifacts (data not shown). 

Using scrambled siRNA-transfected control cells as a reference, the knockdown of most cellular 
kinases was found to be non-cytotoxic within the time frame of this experiment (Fig. 3A and Dataset S1). 
Transfection of siRNAs was considered to be cytotoxic when the viability of cells transfected with a 
target-specific siRNA pool was <85% of the viability of control cells transfected with scrambled siRNAs 
(Fig. 3A). Using this criterion, 222 out of 779 (28.5%) transfections with the specific siRNA pools 
appeared to be toxic to the cells. A minor fraction (50 targets; 6.4%) appeared to be highly detrimental 
(normalized viability value below 75%). To prevent false-positive proviral hits due to a general negative 
effect on cell viability or cell division, we excluded all targets whose knockdown was associated with 
viability measurements below 85%. Such data filtering was not applied for antiviral hits (i.e. hits whose 
knockdown enhanced GFP expression) since siRNA-induced cytotoxicity is expected to inhibit virus 


replication and should therefore not give rise to false-positive antiviral hits. 


Proviral and antiviral proteins and pathways in SARS-CoV-GFP infection. After exclusion of 
cytotoxic siRNA SMARTpools that decreased GFP expression (see above), the remaining 684 targets 
were ranked on the basis of the GFP signal in host factor-depleted SARS-CoV-GFP-infected cells 
compared to control cells (Fig. 3B). Targets were qualified as antiviral or proviral hits if GFP expression 
differed significantly from that in infected control cells transfected with the scrambled siRNA pool (p < 
0.05). Knockdown of the majority of the targets (552 proteins) did not significantly alter GFP reporter 
gene expression (p > 0.05). However, as is not uncommon in this type of screening experiments and 


considering the issue with the origin of the cell line used (see above), we cannot formally exclude that our 
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results were influenced to some extent by insufficient knockdown of certain target genes by the siRNA 
pools in the library. 

Using the criteria outlined above, a total of 90 cellular proteins (19.4% of all targets) were 
identified as antiviral factors, since their depletion significantly increased GFP expression, although for 
most of them less than two-fold. The ten best antiviral hits are depicted in Fig. 4A and the complete 
dataset is provided in Dataset S1. Forty proviral factors were identified and the knockdown of nine of 
those reduced GFP expression by more than two-fold (Fig. 4B; for the complete dataset, see Dataset S1). 
Although, according to the criteria formulated above (p < 0.05), ANGPT4 (214%; p = 0.0555) and PKR 
(210%; p = 0.0884) formally did not qualify as antiviral hits, we have included these proteins in view of 
the exceptionally strong stimulation of GFP expression triggered by their knockdown (Fig. 4A). 
Furthermore, since its knockdown resulted in an almost 3-fold decrease of the GFP signal (35%; p = 
0.0004), DGKE was included as a proviral hit, despite the fact that the viability assay did not rigorously 
exclude cytotoxic effects for this siRNA pool (viability 88%, p = 0.0540). 

The pro- and antiviral hits identified in the siRNA screen were mapped to cellular pathways using 
the IPA software package. Fig. 5 shows the canonical pathways and more general functional categories 
(highlighted in color) in which the proviral (green) and antiviral (red) hits were strongly represented (p < 
0.05). These pathways included apoptosis, cellular immune response, growth factor signaling, cellular 


homeostasis, metabolism of complex lipids, and intracellular and second messenger signaling. 


Evaluation of antiviral hits. An unexpectedly large number of antiviral hits was identified in the primary 
siRNA screen, although for most of them knock-down resulted in a less than 2-fold increase in SARS- 
CoV driven GFP expression. To assess the overall quality of our siRNA screen, and the reliability of the 
identification of antiviral hits in particular, we selected a set of strong and weak antiviral hits for further 
evaluation, namely PKR, ANGPT4, CLK1 (>2-fold increase in GFP signal), CDK6 (1.8-fold increase), 


MAP2K3 (1.6-fold increase), and MAP2K1 (1.2-fold increase). Per target, a deconvoluted set of four 
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individual siRNAs was used for additional knock-down experiments, after which SARS-CoV-GFP 
replication was quantified and the knock-down efficiency at the protein level for each target was evaluated 
by Western blot analysis (Fig. 6). We checked whether there was strong similarity between the target 
sequence of our (human) siRNAs and the corresponding Macaca mulatta sequence, considering the origin 
of the 293/ACE2 cells (see above). This was the case for all siRNAs except for siRNA CDK6 #4, which 
was therefore excluded from further analysis. Using commercially available antisera against the human 
proteins, we were unable to reliably detect endogenous expression of ANGPT4 and CLK1. Therefore, it 
remains uncertain whether ANGPT4 and CLK1 are true antiviral or false-positive hits, as we could not 
determine knockdown levels and correlate these to effects on virus replication (data not shown). 

The mitogen-activated protein kinases (MAPKs) were relatively highly represented among the 
antiviral hits (see Fig. 5 and Supplemental Dataset S1) and therefore we included MAP2K1 and MAP2K3 
in our secondary evaluation. MAP2K1 was a weak antiviral hit in the primary screen, as its depletion led 
to a ~1.2-fold increased GFP expression. In validation experiments with individual siRNAs, we also 
observed a small, but non-significant increase in SARS-CoV-driven GFP expression (Fig. 6A). Western 
blot analysis of siRNA-transfected cells that were infected with wt SARS-CoV revealed poor knockdown 
efficiencies and a clear correlation between the level of MAP2K1 and SARS-CoV N protein expression 
could not be established (Fig. 6B). The siRNA that gave the best knockdown of MAP2K1 (#2) had no 
effect on GFP expression, suggesting that this weak antiviral hit was a false positive in the primary screen. 
Knock-down of MAP2K3 resulted in a >1.6-fold increased SARS-CoV-driven GFP expression in the 
primary screen (see Dataset S1). In cells transfected with a deconvoluted set of individual siRNAs 
targeting the MAP2K3 mRNA we observed a significant increase in GFP expression for 3 out of 4 
siRNAs (Fig. 6C). However, in siRNA-transfected cells that were infected with wt SARS-CoV, only 
introduction of siRNA #3 led to clearly enhanced SARS-CoV N protein expression (Fig. 6D). The other 


siRNAs actually reduced expression of N protein, with an apparent correlation between the remaining 
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percentage of MAP2K3 and N protein levels. Based on these results, MAP2K3 could therefore not be 
confirmed as an antiviral hit. 

We were able to detect expression of cyclin-dependent kinase 6 (CDK6; 1.8-fold increase in GFP 
expression in the primary screen) and found that CDK6 siRNAs #1 and #2 reduced protein levels by at 
least two thirds (Fig. 6F). Transfection with the same siRNAs significantly enhanced SARS-CoV-GFP 
replication (Fig. 6E) and in cells infected with wt SARS-CoV this led to a ~1.5 to 2-fold increase in N 
protein levels (Fig. 6F). These results suggest CDK6 to be a bona fide antiviral hit, as its depletion leads to 
a moderate but significant and reproducible increase in SARS-CoV replication. 

Taken together, the results of our validation experiments suggest that the antiviral hits identified in 
the primary screen should be considered with caution, as several of them may have been false-positives, 


especially those that had a moderate (but significant) effect in the primary screen. 


Validation of PKR as an antiviral factor in SARS-CoV replication. PKR was one of the strongest of 
the 90 antiviral hits that were identified in the primary siRNA library screen. In two independent follow- 
up experiments with re-ordered PKR-specific siRNA SMARTpools, a more than 2-fold increase in GFP 
expression by SARS-CoV-GFP was observed (data not shown), suggesting that PKR is a bona fide 
antiviral hit. PKR is a serine/threonine protein kinase that is activated by double-stranded (ds)RNA, a 
hallmark of RNA virus infection, and the activated form of PKR blocks translation initiation through eIF- 
2a phosphorylation (reviewed in (50)). 

To further validate the antiviral role of PKR in SARS-CoV replication, a deconvoluted set of four 
single PKR-directed siRNAs was used, and transfection of 293/ACE2 cells with three of these siRNAs (# 
2, 3, and 4) significantly increased SARS-CoV-driven GFP expression (Fig. 7A; black bars). Cell viability 
was slightly reduced after transfection with these PKR-directed siRNAs, in particular using siRNA 2 


which caused a 14% reduction in cell viability (Fig. 7A; grey bars). Nevertheless, despite the fact that this 
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siRNA adversely affected cell viability, an increase rather than a decrease of SARS-CoV-driven GFP 
expression was observed. 

Transfection with PKR-specific siRNAs reduced PKR levels in 293/ACE2 cells up to 87% 
compared to control cells, depending on the siRNA used (Fig. 7B). To verify that PKR knockdown 
increased wt SARS-CoV replication, siRNA-transfected 293/ACE2 cells were infected with wt SARS- 
CoV and viral protein expression was analyzed by Western blotting. In line with the effect of PKR siRNA 
#2 on 293/ACE2 cell viability (Fig. 7A), cells transfected with this siRNA contained reduced levels of B- 
actin, which was used as loading control (Fig. 7C; lower panel). Transfection with two of the four 
individual PKR-directed siRNAs (#2 and 3) clearly increased the expression of SARS-CoV N protein 
(Fig. 7C, upper panel), and also led to a ~1-log increase in infectious progeny titers (Fig. 7D). Taken 
together, the increases in GFP signal, N expression and infectious progeny titer correlate well with the 
magnitude of PKR knockdown, which confirms a strong antiviral role for PKR in SARS-CoV-infected 


cells. 


Confirmation of a proviral role for protein kinase C iota in SARS-CoV replication. 

For the evaluation of the proviral hits of the primary screen (see Fig. 4B) four hits were selected that 
caused either a 5-fold reduction (COPB2 and CDK5R2) or a more moderate 2-fold reduction in SARS- 
CoV-driven GFP expression (IHPK1 and PRKC1). We were unable to detect endogenous CDK5R2 and 
IHPK1 by Western blot, and therefore could not validate the proviral role of these two host factors (data 
not shown). The proviral role of PRKC1 and COPB2 could be validated as discussed in the sections below. 
The proviral effect of protein kinase C iota (PRKCt) was validated using the chemical inhibitor sodium 
aurothiomalate (ATM), which blocks the interaction between PRKCt and other PB1-domain-containing 
proteins (51, 52). VeroE6 and 293/ACE2 cells infected with SARS-CoV-GFP were treated with 0.13 to 
20uM ATM, starting 2 h prior to infection. Both in 293/ACE2 (Fig. 8A) and in VeroE6 cells (Fig. 8B), 


SARS-CoV-mediated GFP expression was efficiently inhibited by ATM in a dose-dependent manner with 
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ECs values of 0.58 and 1.06 uM, respectively. No cytotoxicity was observed at the ATM concentrations 


used (Fig. 8). 


COPB2 and proteins of the early secretory pathway are important for SARS-CoV replication. 
COPB2 (or B’-COP) was identified as the strongest proviral hit in our screen, as its knockdown resulted in 
an 82% decrease of GFP expression (Fig. 4B). The coatomer protein complex, of which COPB2 is a 
subunit, contains a total of seven protein subunits (a-, B-, B’-, y-, 6-, ¢-, and C-COP), and drives the 
formation of COPI-coated vesicles, which function in retrograde transport in the early secretory pathway 
(53). To validate its role as a proviral host factor in SARS-CoV replication, COPB2 was depleted by 
transducing 293/ACE2 cells with lentiviruses expressing COPB2 mRNA-specific shRNAs. This reduced 
COPB2 levels by ~70%, compared to control cells transduced with a lentivirus expressing a scrambled 
shRNA (Fig. 9A), a reduction that did not affect cell viability (Fig. 9B). Subsequent infection of COPB2- 
depleted cells with SARS-CoV-GFP resulted in a strong decrease of N protein and GFP expression (Fig. 
9C; left panels). To exclude that the observed effect was an artifact caused by the use of the GFP reporter 
virus, we also analyzed viral protein expression and virus yield in COPB2-depleted cells infected with wt 
SARS-CoV. As for SARS-CoV-GFP, a clear reduction in N protein expression was then also observed in 
COPB2-depleted cells, compared to cells transduced with a lentivirus expressing a scrambled shRNA 
(Fig. 9C; right panels). Titration of culture supernatants from SARS-CoV-GFP-infected cells and wt 
SARS-CoV-infected cells revealed a 2- to 3-log reduction for both viruses upon COPB2 depletion (Fig. 
9D). 

To further substantiate the importance of COPI-coated vesicles for SARS-CoV replication, 
another component of the coatomer protein complex, subunit B1 (COPB1) was depleted by transfection of 
293/ACE2 cells with a COPB1 mRNA-specific siRNA SMARTpool. Depletion of COPB1 resulted in a 
83% reduction of SARS-CoV-driven GFP expression (Fig. 9E). The formation of COPI-coated vesicles is 
mediated through activation of ADP-ribosylation factor 1 (Arfl) by Golgi-specific brefeldin A-resistance 
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430 — guanine nucleotide exchange factor 1 (GBF 1) (54). Therefore, we also analyzed the importance of GBF 1. 
431 GFP reporter gene expression by SARS-CoV-GFP was reduced by 89% in 293/ACE2 cells that had been 
432 depleted for GBF1 (Fig. 9E). GBF1 and COPB1 depletion had no significant effect on cell viability (Fig. 
433 9E). Taken together, these data suggest that COPB2 and COPI-coated vesicles play an essential role in 
434 SARS-CoV replication. 
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Discussion 

In the past decade, functional genomics studies have systematically identified host factors that can 
influence the replication of diverse +RNA viruses (3, 4, 8-10, 39, 55). Here we describe a human kinome- 
wide siRNA screen that aimed to identify factors influencing the entry and replication of SARS-CoV. To 
our knowledge, this is the first report on a systematic functional genomics study of this kind for any 
coronavirus. As kinases are key regulators of many cellular processes, the pro- and antiviral factors 
identified in this study should pinpoint cellular pathways that are important for SARS-CoV replication. 

After we had completed our screen we unfortunately discovered that the 293/ACE2 cells used 
were not of human origin, but must have derived from a non-human primate, probably an Old-world 
monkey closely related to Macaca mulatta. This may have increased the number of false-negative hits, 
due to mismatches between siRNAs designed to target human genes and the sequence of the homologous 
monkey mRNAs. However, the human and Macaca mulatta genome are 94% identical and even a 
nucleotide mismatch in an siRNA would not automatically render it inactive, as it might still silence gene 
expression by blocking translation of the mRNA (56). We therefore concluded that, despite this post- 
screening complication, we should still be able to identify host factors that are relevant for coronavirus 
infection. Moreover, it is important to stress that the use of a non-human primate cell line should not have 
increased the number of false-positive hits. The cell line was highly susceptible to both SARS-CoV and 
EAV infection and could be efficiently transfected with siRNAs. This allowed us to perform siRNA 
screens for host factors involved in the replication of these two distantly-related nidoviruses, and to 
directly compare hits (K.F. Wannee, A.H. de Wilde et a/., manuscript in preparation). 

A recombinant SARS-CoV-GFP reporter virus, in which ORF7a and ORF7b were replaced by the 
GFP gene, was used in our screen in order to conveniently quantify the effect of gene knock-down on 
virus replication. The SARS-CoV ORF7a protein is known to interact with the structural envelope (E), 
membrane (M), and spike (S) proteins (57), and some studies suggest that it may be involved in specific 


virus-host interactions (35, 58, 59). However, the replication efficiency in cell culture of SARS-CoV 
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mutants lacking both ORF7a and ORF7b is unchanged (60, 61), nor were differences in replication 
kinetics, morbidity, and mortality observed in a hamster infection model (61). Although the deletion of the 
two accessory protein genes likely has affected the results of our primary screen only marginally, a wt 
SARS-CoV isolate was used in several of the validation experiments to rule out artifacts caused by the 
lack of expression of ORF7a and 7b. In all cases tested, we did not find major differences between wild- 
type virus and the deletion mutant. 

For SARS-CoV, screening of the kinome-directed library of 779 siRNA SMARTpools resulted in 
the identification of 90 antiviral and 40 proviral proteins. Canonical cellular processes and pathways in 
which these factors were represented strongly include inositol phosphate metabolism, signaling by Rho 
family GTPases, and SAPK/JNK signaling (Fig. 5). Many of the hits could also be mapped to the 
interleukin (IL)-2, -6, -8, and IL-17 signaling pathways, which have previously been implicated in 
controlling coronavirus infection and coronavirus-induced inflammation (reviewed in (2)). For example, 
the SARS-CoV spike (S) protein was shown to induce the expression of the pro-inflammatory cytokine 
IL-8 (62), and IL-6 and IL-8 levels were elevated in the serum of SARS-CoV-infected patients (62, 63). 
Furthermore, mouse hepatitis virus (MHV) and infectious bronchitis virus (IBV) infections were reported 
to upregulate the synthesis of these same cytokines (64, 65). Although our siRNA library screen did not 
target interleukins directly, the identification of (kinase-regulated) interleukin signaling pathways is in line 
with these earlier studies, and emphasizes their importance in SARS-CoV infection. 

A list of host proteins involved in SARS-CoV infection is not more than a good starting point for 
follow-up studies into the role of individual host protein or pathways in the virus replication cycle. 
Previous studies on other viruses, e.g. HIV-1, showed a very limited overlap between hits from 
independent siRNA screens performed in different laboratories (66), highlighting the importance of 
validation and follow-up studies. To judge the overall quality of our siRNA screen, six hits with variable 
impact on SARS-CoV-driven GFP expression were chosen for validation. Four of these, PKR, CDK6, 


COPB2, and PRKC1 (Figs. 6-9) could be confirmed. The weak antiviral hits MAP2K1 and MAP2K3 
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could not be confirmed, as in follow-up experiments knock-down could not be achieved or did not 
convincingly affect SARS-CoV replication, respectively. Interestingly, the diacylglycerol kinase was 
highly represented as hits in the primary screen (6 of 8 targets, Supplemental Dataset S1). Although not 
included in follow-up experiments presented here, a related study that aimed to identify host factors with a 
general effect in nidovirus infection (Fig. 4 and K.F. Wannee, A.H. de Wilde et al., manuscript in 
preparation) confirmed that one of these hits, diacylglycerol kinase epsilon (DGKE), plays a role in the 
SARS-CoV replication cycle. 

MAP2K3, a kinase that acts in the p38 MAPK module, was a moderate hit in our primary screen. 
This MAP kinase signaling pathway is involved in multiple processes like regulation of inflammatory 
responses, cell proliferation, and cell cycle progression (reviewed in (67). This pathway has been 
implicated in the replication of other coronaviruses, but its exact role is still not fully understood. 
Activation of p38 MAPK promotes MHV replication (65) and chemical inhibition of the p38 MAPK 
pathway restricts HCoV-229E replication (68). Overexpression of the SARS-CoV ORF3a (69) and ORF7a 
(35) proteins activates the p38 MAPK signaling pathway, but the role of this pathway in SARS-CoV- 
infected cells remains unclear. Our screen identified several proteins from MAPK signaling cascades, but 
our validation studies suggested that MAP2K3 was a false positive hit. This is supported by the fact that 
the inhibitor SB203580 had no effect on SARS-CoV replication in cell culture (data not shown). This 
compound was previously shown to block HCoV-229E infection in cell culture (68) and to increase the 
survival of SARS-CoV-infected mice through reducing the SARS-CoV-induced inflammatory response 
(70). 

CDK6, a kinase involved in cell cycle progression from G1 to S phase (71), was confirmed as an 
antiviral hit. Depletion of CDK6 results in G1 phase cell cycle arrest. In addition CDK6 is also involved in 
NF-«B signaling and co-regulation of inflammatory genes by binding and activation of the p65 subunit of 
NF-«B (72, 73). Consequently, besides the effect on the cell cycle, depletion of CDK6 might also reduce 


the inflammatory response against virus infection. A recent study by DeDiego ef al. highlighted the 
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relevance of NF-«B-mediated inflammation in SARS-CoV-infected mice (74). Several laboratories 
studied the (antiviral) role of CDK6, and cell cycle progression in general, in coronavirus replication. For 
MHY\V, it has been shown that upon high-MOI inoculation, the virus induces cell cycle arrest in the G0/G1 
phase to promote its replication. In addition, CDK6 is downregulated in MHV-infected 17Cl1 cells (75). 
Similar observations have been made for SARS-CoV, with the N protein limiting cell cycle progression 
by reducing CDK4 and CDK6 kinase activity (76). Overexpression of the SARS-CoV ORF7a protein 
induced cell cycle arrest in the GO/G1 phase, however this was not associated with inhibition of CDK4 
and CDK6 activity (59). In our study, the antiviral role of CDK6 was confirmed and the observed antiviral 
effect is in line with previous studies. 

As pointed out above, our screen yielded a relatively high proportion of antiviral hits, although 
their effect on SARS-CoV replication, while being statistically significant, was generally limited. Based 
on our assessment of some of these moderate hits, at least some of them must have been false-positives. 
Knockdown of PKR had the strongest effect (~2-fold increase in GFP expression) on SARS-CoV 
replication, and this hit could be confirmed independently, as three out of four individual PKR-directed 
siRNAs induced a clear increase in SARS-CoV protein expression and virus yield (Fig. 7C-D). PKR is 
one of four mammalian kinases that can phosphorylate eIF-20 in response to stress signals (the others 
being the PKR-like endoplasmic reticulum kinase (PERK), GCN2, and HRI). Many virus families have 
evolved gene products and strategies to counteract or evade the antiviral action of PKR, highlighting the 
importance of this kinase in the antiviral defense. Previously, it was found that PKR inhibits the 
replication of the coronavirus IBV, as overexpression of a dominant negative kinase-defective PKR 
mutant enhanced IBV replication by almost 2-fold. Furthermore, IBV appeared to (weakly) antagonize the 
antiviral activity of PKR through two independent mechanisms, including a partial block of PKR 
activation (77). Interestingly, MHV-A59 infection in L2 or 17C11 cells did not induce PKR activation and 
the sensitivity of MHV to IFN treatment appeared to be PKR-independent (33, 78, 79). TGEV protein 7 


was shown to counteract PKR activation by binding protein phosphatase 1 to dephosphorylate eIF-20 
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(80), which — in support of our findings - also suggests that PKR is involved in controlling coronavirus 
replication. Krahling et al. showed that PKR was activated in SARS-CoV-infected 293/ACE2 cells, but 
concluded that PKR knockdown did not significantly affect virus replication, despite the fact that a ~1-log 
increase in SARS-CoV titer was observed in their experiments (81). This is in contrast to our PKR knock- 
down experiments, which clearly pointed to an antiviral role for PKR (Fig. 8). In our hands, depletion of 
PKR significantly increased SARS-CoV-driven GFP expression (Fig. 8A), and also enhanced N protein 
expression (Fig. 8C) and the production of infectious progeny (Fig. 8D). This discrepancy cannot be due 
to host cell differences, as the same 293/ACE2 cells were used in both studies (81), and might thus be 
attributed to differences in the experimental set-up, choice of controls, or normalization and interpretation 
of the data. 

In line with the findings for PKR, reducing the expression of PERK (or EIF2AK3), one of the 
other kinases known to phosphorylate eIF-2a, resulted in an increase of SARS-CoV-GFP reporter gene 
expression with 57% (p < 0.01; Dataset S1). The unfolded protein response - i.e. the detection of 
misfolded proteins within the ER lumen - activates PERK, which in turn phosphorylates eIF2a, and 
ultimately triggers apoptosis. The relatively strong antiviral effect of PERK observed in this study is in 
line with previous studies suggesting that the phosphorylation of eIF-2a in SARS-CoV-infected cells is 
mediated by PERK activation (81). Our findings support the hypothesis that upon SARS-CoV infection 
the unfolded protein response is activated as an antiviral strategy. In fact, countering SARS-CoV infection 
may involve multiple cellular responses that induce apoptosis, including the activation of PKR and PERK, 
which could also explain the identification of several other hits involved in apoptosis, like those from the 
SAPK/JNK pathway. 

Among the proviral hits, PRKCt had a relatively moderate effect in the primary screen (Fig. 4B), 
but its proviral role was validated using the inhibitor ATM (Fig. 8). Members of the protein kinase C 
family are serine/threonine protein kinases and involved in several signalling pathways that regulate e.g. 


cell proliferation, cell cycle progression, and cell survival (reviewed in (82)). Interestingly, PRKCt can be 
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activated by phospho-inositol lipids involved in microtubule dynamics within the early secretory pathway 
(83). PRKCt contains an N-terminal PB1 domain that ensures the signalling specificity (reviewed in (84)) 
and ATM affects the interaction of PRKC1 with other PB1 domain-containing proteins like Par6é, MEK5, 
and p62 (51, 52, 84). Therefore, blocking the PRKC1 PB1 domain could decrease MEKS5 (85) and NF-kB 
signalling (via p62), and affect cell polarity (via Par6) (reviewed in (84)). In addition, PRKC1 plays an 
essential role in COP-I vesicle formation, since the GTPase Rab2 binds PRKC1 and ultimately promotes 
recruitment of B-COP to pre-Golgi membrane structures for the formation of early secretory vesicles (83, 
86, 87). As discussed below, COPB2 and the early secretory pathway play a crucial role in SARS-CoV 
replication, and in this manner PRKC1 may affect SARS-CoV replication as well. Tisdale et al. have 
shown that PRKCt kinase activity is essential for the generation of retrograde-transport vesicles (87). 
However, the role of the PRKCt PB1 domain, the main target of the drug ATM, in COP-I vesicle 
formation was not directly investigated in this study. Nonetheless, this hypothesis is substantiated by the 
ATM concentration that blocked SARS-CoV replication. Our ECs values are similar to the reported ICs 
(~1 pM) for inhibition of the PRKC1 PB-mediated interactions (51), while the ICs for inhibition of the 
kinase activity was ~100-fold higher (88). The exact mechanism by which PRKCt influences the SARS- 
CoV replicative cycle remains an interesting topic for future research, and it could even be an interesting 
target for the development of host-directed antiviral therapy for pathogenic coronaviruses. 

Coronavirus replication is associated with a cytoplasmic reticulovesicular network of modified 
ER, including double-membrane vesicles and convoluted membranes (28). Despite the in-depth 
characterization of their ultrastructure, the biogenesis of these membrane structures and the cellular factors 
involved have remained largely uncharacterized. For example, the membrane source of these coronavirus- 
induced replication structures is still controversial, with advanced EM analyses showing the RVN to be 
derived from and continuous with the ER (28, 44, 89) and other studies implicating the autophagy 
pathway (90) or EDEMosomes (91) as the primary membrane donor. Our earlier work already suggested 


that the integrity of the early secretory pathway is important for efficient SARS-CoV replication, as 
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brefeldin A treatment of SARS-CoV-infected cells significantly reduced replication as well as the 
accumulation of virus-induced membrane structures (89). In line with these findings, COPI-coated 
vesicles were implicated in the biogenesis of MHV replication structures (92, 93) and SARS-CoV nsp3 
was shown to interact with three COPI subunits (94). In none of these previous SARS-CoV and MHV 
studies a complete block of virus replication could be achieved, neither by reducing COPI vesicle 
formation by depletion of one of the coatomer subunits, nor by treatment with brefeldin A. These results 
may in part be explained by incomplete knockdown or the presence of residual COPI vesicles (complete 
knockdown is probably not possible due to its detrimental effect on intracellular trafficking and cell 
viability). Although our present study clearly demonstrates the importance of COPI-vesicles in SARS- 
CoV replication (Fig. 9), their role in the formation or function of the SARS-CoV-induced RVN remains 
elusive. The importance of COPI-coated vesicles is further supported by their essential role in the 
replication of many other RNA viruses, such as poliovirus (95, 96) and other picornaviruses (40, 97-100). 
In conclusion, our kinome-wide siRNA screen has identified several cellular proteins and 
pathways that influence SARS-CoV replication and possibly coronavirus infections in general. Our data 
thus provide a starting point for further validation and in-depth mechanistic studies which should enhance 


our understanding of the complex interplay between coronaviruses and their host. 
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Figure legends: 


Fig. 1. Viability and susceptibility to SARS-CoV infection of siRNA-transfected 293/ACE2 cells. (A) 
293/ACE2 cells were infected with SARS-CoV-GFP (MOI 10) and at 24, 28, and 30 h p.i. cells were 
fixed and GFP fluorescence was measured. (B) 293/ACE2 cells were transfected with siRNAs targeting 
GAPDH mRNA or a scrambled control siRNA (Scr). At 48 h p.t., cells were infected with SARS-CoV- 
GFP (MOI 10) and 24 h later cells were fixed and GFP expression was measured (black bars). Cell 
viability (dark grey bars) was analyzed at 48 h after siRNA transfection and knockdown of GADPH 
expression was monitored by measuring enzymatic activity (light grey bars). All values were normalized 


to those obtained with non-transfected control cells (100%). 


Fig. 2. Design of siRNA library screening procedure. See text for details. 


Fig. 3. Results of the siRNA screens for host factors influencing SARS-CoV replication. (A) Viability 
assays were done at 48 h p.t. and data were normalized to the measurements for control cells transfected 
with scrambled siRNA (100%). Data were binned into 4 viability categories as indicated below the x-as 
and the number in each bar is the absolute number of siRNA targets within that category. The fraction of 
the total (779) number of targets in each category is indicated above each bar. For each siRNA pool in the 
library, the viability data are the average of nine measurements, resulting from three independent library 
screens. (B) The plot shows the distribution of the log.-transformed values of GFP reporter gene 
expression by SARS-CoV-GFP in siRNA-transfected cells, normalized to the GFP signal of infected 
control cells that were transfected with scrambled siRNA. Targets were ranked based on the magnitude of 
the effect of their knockdown on SARS-CoV replication. Targets were considered to have a robust 
antiviral effect when their knockdown increased reporter gene expression to at least 150% (red area above 


x-axis). Proviral hits, whose knockdown induced an at least 2-fold reduction in GFP expression, are 
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depicted in the green area below the x-axis. Proviral targets whose knockdown reduced cell viability to 
below 85% were excluded (see main text), leaving a total of 684 targets included in the final analysis. The 


plot represents the average of three library screens (each done in triplicate). 


Fig. 4. Heat-maps of pro- and antiviral hits identified in this study. (A) List of the ten most prominent 
antiviral (A) and proviral hits (B). For each target, the p-value, accession number, and gene name are 
shown. Each data point represents the result of a single library screen and is the average of the 3 replicates 


that were done in each screen. The full hit-lists can be found in Dataset S1. 


Fig. 5. Cellular pathways influencing SARS-CoV-GFP replication. Graphical representation of the 
canonical pathways (white ellipses) identified in the siRNA library screen for cellular factors affecting 
SARS-CoV replication. All proviral (green) and antiviral hits (red) (for a complete list, see Dataset $1) for 
which depletion significantly altered SARS-CoV-GFP replication were used to identify cellular pathways 
by IPA in which the hits were clearly overrepresented, and only those pathways and the hits represented in 
them are shown here. The hits are represented by nodes with lines linking them to one or more canonical 
pathways. The color intensity of the nodes indicates the strength of the pro- or antiviral effect (log, ratio of 
GFP expression normalized to infected control cells), e.g. factors with a stronger antiviral effect have a 
more intense red color. The identified canonical pathways were clustered into more general categories that 


are indicated by text boxes in the colored background shading. 


Fig. 6. Evaluation of the antiviral hits CDK6, MAP2K1, and MAP2K3. 293/ACE2 cells were 
transfected with four individual siRNAs targeting, MAP2K1 (A, B) or MAP2K3 (C, D), or three CDK6- 
specific siRNAs (E, F). A non-targeting scrambled siRNA was used as a control. At 48 h p.t. cells were 
infected with SARS-CoV-GFP at an MOI of 10 (A, C, E), fixed 24 h later, and GFP fluorescence (black 


bars) was quantified and normalized to the value measured in infected, scrambled siRNA-transfected cells 
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(100%). The effect of siRNA transfection on cell viability was analyzed in parallel (grey bars) and values 


were normalized to those of scrambled siRNA-transfected control cells (100%). Each experiment was 


repeated at least three times (average + SD). In parallel, siRNA-transfected cells were infected with wt 
SARS-CoV (MOI 5), and at 8 h p.i., SARS-CoV N expression was monitored by Western blotting (B, D, 
F). TfR was used as a loading control. Knockdown levels of the host proteins were analyzed by Western 
blot (B, D, F). The amount of SARS-CoV N protein and remaining quantity of host protein compared to 
scrambled siRNA-transfected cells (100%) is shown below each lane. All experiments were repeated at 


least twice. 


Fig. 7. Validation of PKR as an antiviral factor in SARS-CoV replication. 293/ACE2 cells were 
transfected with four individual siRNAs targeting PKR or a scrambled control siRNA. (A) At 48 h p.t. 
cells were infected with SARS-CoV-GFP (MOI 10), fixed 24 later, and GFP fluorescence (black bars) was 
quantified and normalized to the value measured in infected, scrambled siRNA-transfected cells (100%). 


The effect of siRNA transfection on cell viability was analyzed in parallel (grey bars) and values were 


normalized to those of scrambled siRNA-transfected control cells (100%). Average + SD is given (""; p- 
value < 0.001). (B) Knockdown of PKR expression at 48 h p.t. was monitored by Western blotting and the 
percentage of PKR remaining compared to scrambled siRNA-transfected cells is shown below each lane. 
TfR was used as loading control. (C) Cells transfected with PKR-specific siRNAs and control cells were 
infected with SARS-CoV (MOI 0.01) and 24 h later these cells were lysed to assess SARS-CoV N levels 
by Western blotting (shown below the panels as percentage of control), using B-actin as loading control. 
(D) Virus titers in the 24-h p.i. culture supernatants of wt SARS-CoV-infected cells (MOI of 0.01) 


transfected with PKR-specific or scrambled siRNA. All experiments were repeated at least twice. 


Fig. 8. Validation of PRKC. as proviral host factor. (A) 293/ACE2 or (B) VeroE6 cells in 96-well 


plates were infected with SARS-CoV-GFP (MOI 10). Treatment with 0.13-20 uM sodium aurothiomalate 
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(ATM) was started 2 h prior to infection and cells were fixed at 18 h p.i. (VeroE6) or 24 h p.i. 
(293/ACE2). GFP reporter gene expression (black bars) was measured and normalized to the signal in 
untreated control cells (100 %). The grey lines show the effect of ATM on the viability of mock-infected 
cells, normalized to the viability of solvent-treated control cells. Graphs show the results (average and SD) 


of a representative experiment performed in quadruplicate. Both experiments were repeated at least twice. 


Fig. 9. Proteins of the early secretory pathway are important for SARS-CoV replication. (A) 
293/ACE2 cells were transduced with lentiviruses expressing a COPB2 mRNA-specific or a scrambled 
shRNA. Knockdown of COPB2 expression at 48 h p.t. was monitored by Western blotting with a COPB2- 
specific antiserum and cyclophilin B (CypB) was used as loading control. (B) Viability of COPB2- 
depleted 293/ACE2 cells was analyzed at 48 h after transduction (% of control cells transduced with 
lentiviruses expressing a scrambled shRNA). (C, D) COPB2-depleted and control cells were infected with 
either SARS-CoV-GFP or wt SARS-CoV (MOI of 0.01). (C) SARS-CoV protein expression at 32 h p.i. 
(SARS-CoV-GFP) or 24 h p.i. (wt SARS-CoV) was analyzed by Western blotting with N-specific and 
GFP-specific antisera, using the TfR protein as loading control. SARS-CoV N protein expression was 
quantified and normalized to that inscrambled siRNA-transfected cells (100%) as indicated under each 
lane. (D) SARS-CoV-GFP (black bars) and wt SARS-CoV (grey bars) progeny titers in the culture 
supernatants of control or COPB2-depleted cells at 32 h p.i. (SARS-CoV-GFP) or 24 h p.i. (wt SARS- 
CoV). (E) Normalized GFP expression by SARS-CoV-GFP in 293/ACE2 cells transfected with siRNA 
SMARTpools targeting COPB1, GBF1, or a scrambled control siRNA. Cells were infected 48 h p.t. at an 
MOI of 10 and 24 h later GFP fluorescence was quantified and normalized to that in infected cells 
transfected with a scrambled siRNA. GFP fluorescence data is the average of three independent 


experiments. 
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A GFP expression 
(fold increase) 
GeneID 1 2 3. average p-value Accession 
ANGPT4 ; : : : NM_015985 
PKR ! : ; 2.1 0.0884* NM_002759 
CLK1 : : : 2.0 0.0077 NM_004071 
MAP2K6 : : : 1.9 0.0014 NM_002758 
CSNK1G1 ; : E 1.9 0.0067 NM_022048 
EPHA3 : : : 1.8 0.0176 NM_005233 
CDK6 : , : 1.8 0.0140 NM_001259 
AURKB : ! : 1.8 0.0016 NM_004217 
GCK : : : 1.8 0.0119 NM_000162 
DGKD : : : 1.8 0.0135 NM_003648 
B GFP expression 
(fold increase) 
GeneID 1 2 3 average p-value Accession 
COPB2 0.2 0.0143 NM_004766 
CDK5R2 0.2 0.0251 NM_003936 
PFTK1 0.3 0.0087 NM_012395 
ABI1 0.3 0.0055 NM_005470 
DGKE 0.4 0.0004" NM_003647 
NME2 0.4 0.0027 NM_002512 
AZU1 0.4 0.0075 NM_001700 
IHPK1 05 0.0405 NM_153273 
PSKH1 0.5 0.0094 NM_006742 
PRKCI 0.5 0.0157 NM_002740 
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* Note: p>0.05, but PKR and ANGPT4 were also included as antiviral hits based on their large effect on reporter gene expression 
T Note: siRNAs have a cytotoxic effect (88% viability, p = 0.0540) 
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Figure 6 
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